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Statistical properties of nonlinear systems revealed in an unexpected way already 
in Fermi-Pasta-Ulam paper, are responsible for the appearance of substructures. 
We represent a method within S-tree diagram formalism for investigation of the 
substructures of N-body gravitating systems in a defined time interval. This ap- 
proach enables one to consider the link between substructures of the system and 
its relative instability and is efficient for the study of the substructures in clusters 
of galaxies. 

1 Introduction 

The study of hierarchical distribution of galaxies is one important issues of 
observational cosmology since the properties of substructure of galaxy clusters 
have to carry information on the mechanisms of their formation and evolution. 

Various methods are developed for the study of the substructures in 
galaxy clusters, e.g. two and higher-point correlation functions, Li statis- 
tics, topological measures, wavelets, etc which are using various descriptors 
and assumptions Llffla. S-tree method □ is using both the positional and kine- 
matical information about the system in self-consistent way, by means of the 
consideration of the geometrical properties of the phase space of A^-body sys- 
tem and the definition of the criterion of boundness between particles. That 
method enabled to reveal the substructures in various galaxy clusters, i.e. 
Perseus, Virgo n Coma, Abell ENACS clusters, etc HCteI. The existing versions 
of the S-tree U'H are dealing with the substructures at the fixed time. 

In this paper we will study the behavior of substructure within certain 
time interval i.e. to try to trace the future and past evolution of the subsys- 
tems within that interval. 

The approach is based on the generalized S'-scheme, and therefore, first, 
we will briefly review the main steps of the approach. 

2 N-body bounded systems: 5-diagram method 

The key idea of the method, as already mentioned above, is the introduction 
of the concept of the degree of boundness p. 
Consider a set of N points 
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X = {xi, . . .,x N }, 

the function P 

P : X x X -> R + and p G R+. 

The basic definitions are as follows. 
Definition 1 

We say that Vi 6 X and Vy 6 X are p-bounded, if P(ic, y) > p. 
Definition 2 

We say that U C X(U ^ 0) is a p-bounded subgroup, if: 

1. Vx G {7 and Vy £ U P(x, y) < p; 

2. \/x E U and \fy G U 3x = Xi x , X{ 2 , . . . , Xi h = y, that P(xi l , Xi l+1 ) > 
p; Vl = l,...,k-1. 

Definition 3 

We say that U\, . . . , Ud is the distribution of the set X via p-bounded 
groups, if: 

i- Uti Ui = x- 

2. i^j (i,j = l,...,d)^U i r\U j = <t>; 

3. Ui(i = 1, . . . , d) is a p-bounded group. 

The function P Different physical quantities can be assigned to the function 
P as listed in such as the distances, the energy, the force, the perturbations 
of potential, momentum, and so on. For example, the choice of the mu- 
tual distance of the particles in order to define the subgrouping is equivalent 
to the corresponding correlation functions (spatial or angular, related with 
each other via the Limber equation), which is, obviously, a rather incomplete 
characteristic for that aim. It can be shown that, at least for astrophysical 
problems, among the most informative ones is the Riemannian curvature of 
the configuration space El, which determines the behavior of close geodesies, 
as known from basic courses on classical mechanics E3. 

So, by this, S-tree algorithm we obtain the substructure of the N-body 
system, i.e. the degree of boundness of each subgroup of particles, for any 
given function P and Vp. This found distribution will satisfy Definitions 1, 
2 and 3. The final result can be represented either through tables or by 
tree-diagrams, graphs (S-tree). 



bekar: submitted to World Scientific on February 2, 2008 2 



The next important aspect is that the method allows to use simultane- 
ously various boiindness criteria between the particles of the system for fixed 
moment of time u. 

Consider the iV-body system: X — {xi, x%, . . . , xn} and Pi, ... , Pjj func- 
tions. Consider also Di, . . . , Dh matrices, where D a = (dfj), a = 1, . . . , H; 
i,j = 1, . . . , N; d?i = 0. If i y£ j = P a (xi,Xj). Constructing matrix D 
in the following way: 

r> = (4), i,j = i,...,N, 

where 

dij = {dlpdl, ...,df } ), i,j = l,...,N. 

D matrix contains the whole information about matrices D\, . . . ,Dh- For 

any function P a (a = 1, . . . , H) exists its own boundness vector 

Pa = [pi j ■ ■ ■ ! Pq a )- We construct for the D matrix a vector of boundness 

Pd = (p?,---,Pq ) 

and for k = 1, . . . ,Q D ] Pk= (Pki> • ■ ■ >/°fc^)) wnere k a = l,...,Q a and a = 
1, . . . , H. Vector pn is constructed by all possible combinations of components 
of pf? vector. For any function P a , (a = 1, . ■ ■ , H) we introduce "the degree 
of influence" W a G . The next steps are realized for fixed value of p® . We 
transit from the matrix D to the matrix D u in a following way 

D u = (uij); i,j = l,...,N, 

where My = (u\- , . . . , u 1 ?) and ufj, a = 1, . . . , H is defined as 

„a = {W a , if dg > PL 

The following D v matrix we construct as 

H 

D v = (Vij), i, j = 1, • ■ ■ ,N, v i:j = ^2 u fj- 

a=l 

Defining p = (p 1; . . . ,ps), the vector of boundness for the D v matrix and, 
using the S'-method for D v and for any given component of p vector, we 
obtain the distribution of the initial system on the bounded subgroups. Note 
finally, that according to the ^-generalized scheme with fixed values of p® 
and p;, we obtain the distribution of system on p-bounded subgroups, where 
l = l,...,S; p={pk\m)- 

Consider a dynamical time interval or including z dynamical interval 
[ta,fb] of an initial system Eil. Divide this interval on H — 1 equal parts. 
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For any t a ; a = 1, . . . , H it is possible to find all sets of (Li(t a ), Vi(t a )), and 
where i = 1, . . . , N; Li,Vi - coordinates and velocities of i-th particle. Con- 
sider the P function at different moments of time. Represent the function 



Pi, . . . , P H as Pi = P(h) = P(t a ), P 2 = P(t 2 ), ...,Ph= P{t H ) = P(t b ). 



As in the previous section, we will obtain D and D u matrices. 

For W a — 1, a = 1, ... H, W a — 1. For the matrix D v and /i G R + using 
the .S-tree method we obtain the splitting of the system into on /it-boundcd 
subgroups within the about time interval [i a ,£&]- 

As a prefered floor of distribution we chose the larger components of 
the vector /2, because /2-connection is a quantitative estimation of the sequence 
of numbers, where 



After splitting the system into bounded groups by the suggested model, $/7V 
value may present a define physical interest, jszhere $ - is numebr of those 



Consider two different systems Xi,X 2 with Ni and N 2 particles. For 
both systems a distribution is analyzed, and [^i^b] intervals have to 

be considered in a self- consistent way. The difference of $i/iVi and & 2 /N 2 
can be considered as a criterion of relative instability of these systems. Note, 
that the presented approach allows to link directly the substructure properties 
and the relative instability of the N-body systems. 

In order to apply this scheme for physical systems the correspondence 
between the chosen time intervals [t^, t\] and \t 2 al tf\ have to be defined taking 
into account the formulation of the criterion of relative instability for systems 
X\ and X 2 . 

3 Conclusion 

The study of the substructure of nonlinear N-body systems is an interesting 
theoretical problem, as well has a direct link with the astrophysical systems. 
Particularly, the S-tree approach had proved its practical efficiency in the 
study of the substructure of clusters of galaxies. Particularly, the subgroups 
of the Coma cluster detected by S-tree □ have been later studied observa- 
tionally and certain predicted trends in the properties of galaxies of one of 
subgroups (subgroup 2) have found their apparent confirmation 113. The fur- 
ther development of the method therefore, seems of particular interest. 

Here we proposed a generalization of the approach enabling the investiga- 
tion of evolution of the substructuring properties of a nonlinear system within 
a defined time interval. Application of this scheme to the clusters of galaxies 
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< = 1, 1 V. a=l,...,H 




id to a valuable insight on the evolution and maybe on the origin of the 
substructures. 

We are grateful to Prof.V.G.Gurzadyan for discussions and advice. 
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